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ABSTRACT 
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o 

D ■ The detection by the HESS atmospheric Cherenkov telescope of fifteen new 



sources from the Galactic plane makes it possible to estimate the contribution 
of unresolved point sources like those detected by HESS to the diffuse Galactic 
emission measured by EGRET and recently at higher energies by the Milagro 
Collaboration. Assuming that HESS sources have all the same intrinsic lumi- 
nosity, the contribution of this new source population can account for most of 
the Milagro 7-ray emission at TeV energies and between 10 and 20 per cent of 
EGRET diffuse Galactic 7-ray emission for energies bigger than 10 GeV. Also, 



by combining the HESS and the Milagro results, constraints can be put on the 
distribution and the luminosities of gamma ray emitters in the Galaxy. 

Q-f 

q ■ Subject headings: gamma rays: theory 



1. Introduction 

The Galactic diffuse 7-ray emission is believed to be mostly produced in interactions 
of cosmic rays with the matter and the radiation fields in the Galaxy, the main production 
mechanisms being electron non-thermal Bremsstrahlung, Inverse Compton scatterings off the 
radiation fields and pion decay processes in inelastic collisions of nuclei and matter. Although 
the standard production mechanisms of 7-rays (Bertsch et al. 1993) explain generally well 
the spatial and energy distribution of the emission below 1 GeV, the model does not match 
EGRET observations of the 7-ray sky above 1 GeV. (Hunter et al. 1997) The discrepancy 
between the model and the observations is more pronounced towards the direction of the 
inner Galaxy, where it amounts to 60 per cent of the observed emission. Many possible 
explanations have been proposed to account for the GeV excess. Aharonian & Atoyan 
(2001) have proposed an hardening of the cosmic ray proton spectrum. The average proton 
spectrum in the Galaxy would have to be harder than locally measured. Strong et al. (1999, 
2004) have suggested a renormalization of the electron spectrum in order to fit EGRET data. 
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Kamae et al. (2006) have argued that diffractive processes make the 7-ray spectrum harder 
than the incident proton spectrum by 0.05 in power-law index, that the scaling violation 
produces more 7r°, and, that the inelastic p — p cross section has a logarithmic increase with 
the incident proton energy, so that the combination of these effects could explain about 20 
per cent of the " GeV Excess" . de Boer et al. (2005) claimed that the features of the EGRET 
excess are consistent with the hypothesis of a signal from Dark Matter Annihilation. Recently 
Strong (2006) has argued that the excess in the diffuse emission measured by EGRET and 
COMPTEL could be explained by a combination of Galactic source populations too dim to 
be detected by EGRET in addition to the standard production mechanisms. 

In any case a major difficulty when studying the diffuse emission is to disantangle the 
truly diffuse emission from that produced by unresolved sources. 

At TeV energies, Milagro has observed diffuse emission from a large region of the Galac- 
tic plane from Galactic longitude / of 40° < / < 100° (Atkins et al. 2005). The different 
theories to explain the GeV excess make very different predictions about the TeV flux. Re- 
cently Prodanovic et al. (2006) have argued that no strong signal of pion decay is seen in the 
7-ray spectrum. The pion decay mechanism would not be able to explain the excess above 
1 GeV. Therefore, Prodanovic et al. (2006) claimed that the Milagro measurements of the 
diffuse flux at higher energies revealed a new excess at TeV energies and concluded that a 
different mechanism is required to explain the two excesses, either dark matter decay or con- 
tribution from unresolved EGRET sources. An extrapolation of the EGRET sources within 
the range in Galactic longitude of the Milagro observation overestimates the diffuse TeV flux 
measured. However, there could be a population of sources undetectable by EGRET, but 
contributing to both the GeV and TeV excess diffuse emission. Recent observations by the 
TeV observatory, HESS, in fact point to a new class of hard spectrum TeV sources. HESS is 
an array of telescopes using the imaging atmospheric Cherenkov technique to detect 7-rays 
above 200 GeV. HESS scans the inner part of the Galaxy in the region between —30° and 
30° degrees in Galactic longitude (—30° < I < 30°) and between —3° and 3° degrees in 
Galactic latitude (—3° < b < 3°) with an angular resolution better than 0.1 degrees. Its 
average sensitivity at the energy threshold of 200 GeV is 2 per cent of the Crab flux within 
±2° degrees in Galactic latitude and decreases outside this range. (The flux of the Crab 
nebula is about 2.4 x lO^dn^s^ 1 for E > 200 GeV (Aharonian et al. 2004)). HESS has 
previously reported the detection of 8 unknown sources close to the Galactic Plane. (Aha- 
ronian et al. 2005) Recently the total number of discovered sources has been updated to 
15. (Aharonian et al. 2006) Three sources have also been detected that are identified with 
the Galactic Center, the SNR RX J1713.7-3946 and GO.9+01. All HESS sources are located 
close to the Galactic Plane, within ±1 degrees of latitude. (Aharonian et al. 2006) The fluxes 
above 200 GeV from all sources are fit to power laws, 4tj = F E~ r , whose photon indices T 
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range from 1.8 to 2.8, with a relatively hard mean index of 2.32 and RMS of 0.2. 

In the inner Galaxy region observed by HESS the existence of 91 SNRs and 389 pulsars 
was already known at lower energies (Green 2004; Manchester et al. 2005). Many of these 
SNRs and pulsars within the inner Galaxy might emit VHE 7 rays but only a few 7 ray 
sources were previously known. Although the HESS Collaboration indicates a clear positional 
coincidence of its sources with a known SNR or pulsar only in a limited number of cases, the 
distribution of Galactic latitude of the eighteen VHE 7 ray sources detected by HESS agrees 
quite well with the distributions of all SNRs catalogued by Green (2004) and of all pulsars 
catalogued by Manchester et al. (2005). Assuming, for example, that SNRs are a single class 
of counterparts with isotropic luminosity 1.95 x 10 34 ergs _1 to the new HESS sources and 
taking for them a simple radiative model, Aharonian et al. (2006) found that the location of 
these sources favours a scale height of less than 100 pc, consistent with the hypothesis that 
these sources are either SNRs or pulsars in a massive star forming region. 

The HESS detections of high energy 7 rays from fifteen new sources has improved sig- 
nificantly the knowledge of both the spatial distribution and the spectra and fluxes of VHE 
7-ray galactic sources. These HESS results make it possible to estimate with unprecedented 
precision the contribution of unresolved point sources to the Galactic diffuse emission mea- 
sured up to 100 GeV by EGRET and recently extended to TeV energies by Milagro. (Atkins 
et al. 2005) Based on HESS results, knowing the sensitivity and the field of view of an ex- 
periment, we estimate the number of expected sources and their expected VHE 7-ray flux. 
Finally we evaluate the contribution of unresolved point sources to the diffuse 7-ray emission 
for EGRET and Milagro. 



2. Method 

Under the simplistic assumption that all sources have the same intrinsic energy lumi- 
nosity L for E > E th , the integral photon flux F from a source depends upon the distance 
D of the detector as 

r-i £;,4iD ! k ' 

where E t h is the detector threshold energy and T is the spectral index if the 7-ray emission 
is a power law. Inverting Eq. (1) 

where we have defined k as 
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If we assume that all sources are located in the Galactic plane and their surface density, 
which does not have to be constant, is a sources , then by differentiating Eq.(2) 

dD = ^ F -^ dF (4) 

we obtain the number of sources per interval of Galactic longitude dl and per interval flux 
dF 

1 k 

sources ® sources dl D dD — <J sources dF dl . (5) 

2 r z 

As first case we assume that cr sources is constant. Then integrating Eq.(5) over all source 
fluxes above HESS sensitivity of about / m j„ = 4.8 x 1CT 12 photons s -1 cm" 2 and over the 
range in longitude / surveyed by HESS the number of sources detected by HESS fixes the 
density of sources in the Galactic plane 



N 

1 ' sources 



dF / Jl^p± = 18, (6) 

.jmin ?T 



where the maximum flux detected by HESS was f max = 6.0 x 10 11 photons s 1 cm 2 . 
For a given source luminosity L from Eq.(6) 



5-7 x 10- 10 . _ 2 _ 

^Lo! k P C • ( 7 ) 
n T-l E th 4tt 



For instance, assuming that the emission follows a power law with spectral index Y = 2.32, 
the average spectral index of HESS sources, if L = 10 33 erg/s, o sources = 5.3 kpc~ 2 , for 
L = 10 34 erg/s, Usources = 0.5 kpc~ 2 . The average distance is 5.6 kpc for this luminosity. 
This range of surface density is intermediate between the local density of SNRs and pulsars. 

The region in the Galaxy where point sources are likely located extends from D m i n = 0.3 
kpc up to D max = 30 kpc (Swordy 2003). Assuming that the flux from each source depends 
only upon the distance from the detector F = F(D) and not on the type of source, HESS 
total flux per steradian from resolved and unresolved sources is 

F(E > 200 GeV, \l\ < 30) = / F(D) dN sources 

= / / KCJ ^es dlDdD 

J J->mm J I — — s 



8.0 x 10 9 photons cm £ s 1 sr^ 



(8) 
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In Eq. 8 the total integral flux is already divided by the HESS field of view. This flux is 
independent of L . 

The surface density of 7-ray sources is not likely to be constant. The distribution is 
probably similar to that of 7-ray candidate sources such as pulsars or supernova remnants 
which are more prominent in the inner galaxy. However, the ratio of 7 emitting pulsars 
versus SNRs is unknown. The two extreme cases of all HESS sources being either pulsars 
or SNRs are considered. <J sources in Eq.(5) will be either the surface density of SNRs or the 
surface density of pulsars in the Galactic Plane, as observed in the radio wavelengths. We 
do not consider the latitude distributions of pulsars and SNRs because most point sources 
are located close to the plane of the Galaxy, within a region (—2° < b < 2°) which extends 
from D min = 0.3 kpc up to D max = 30 kpc (Swordy 2003). In terms of the heliocentric 
distance D and the longitude angle / the galactocentric distance, r, is 

r = \J rl + D 2 — 2r Q D cos (I) (9) 

with r = 8.5 kpc, which is the distance from the Earth to the Galactic center. The pulsar 
surface density o- pu isar(j'), plotted in Fig.l, is fitted by the following shifted Gamma function 
(Yusifov k Kiigiik 2004; Lorimer et al. 2004) 

X a [ b( x ~ x o )] 

Vpulsar(r) = (J 0p ( — ) e l -0 (10) 

where x = r + n and x = r + n, a 0p = 37.6 ± 1.9 kpc~ 2 , a = 1.64 ± 0.11, b = AM ± 0.24 
and r 1 = 0.55 ± 0.10 kpc. The SNR surface density, plotted in Fig.l, is (Green 2004; Case 
& Bhattacharya 1998) 

f a 0snr sin(^ + e )e-P r ior r < 16.8 

<?SNR\r) = < (11) 
[ for r > 16.8 

with a 0snr = 1.96 ± 1.38 kpc~ 2 , r 2 = 17.2 ± 1.9 kpc, 6 = 0.08 ± 0.33 and f3 = 0.13 ± 0.08. 
Integrating Eq.(5) over all pulsar fluxes above HESS sensitivity, f m in, and below f max , the 
greatest flux from a source detected by HESS, and assuming a given source luminosity L , 
the number of sources resolved by HESS in the region (—30° < / < 30°) fixes the maximum 
fraction of pulsars that are 7-ray emitters 

Nhess = J dDD J dla pulsar {l,D) = 18 

^ rfmax f l =J 

= 2 dF F^ n dla pulsar (l,F) = 18. 

^ f min "6 

(12) 
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Galactocentric Radius (kpc) 

Fig. 1. — Pulsar and SNR surface distributions versus galactocentric radius r from Yusifov & 
Kiigiik (2004) and Lorimer et al. (2004) and from Green (2004) and Case & Bhattacharya 
(1998), respectively. 



If HESS population of sources consists entirely of pulsars, the total flux from resolved and 
unresolved sources divided over HESS field of view is 



W(B > 200^,1.1 < 30) = /F^diW 

rD max r l =% 

= dDD— dla pulsar (l,D) 



= 1.0 x 10 8 photons cm 2 s 1 sr 1 . 

(13) 

If the all radio loud pulsars emit 7-rays then the luminosity predicted in Eq. (12) is 1.2 x 10 32 
ergs/s and the average distance would be 0.5 kpc. Though the ratio between 7-ray loud versus 
7-ray quiet pulsars is uncertain. If we suppose that only five per cent of the radio loud pulsar 
emit 7-rays, then the 



pulsar 



(E > 200 GeV, \l\ < 30) = 7.0 x 10~ 9 photons cm" 2 s" 1 sr" 1 , (14) 



the source luminosity predicted would be 1.7 x 10 33 ergs/s and the average distance for HESS 
sources 1.8 kpc. If 1 per cent of radio loud pulsars emit TeV 7s then 

F P uisar(E > 200 GeV, \l\ < 30) = 5.4 x 10~ 9 photons cm" 2 s" 1 sr" 1 . (15) 

the source luminosity predicted would be 6.5 x 10 33 ergs/s and the average distance for HESS 
sources 3.5 kpc. 



-7- 



If instead all HESS sources belong to the SNR population the total flux would be 
F snr {E > 200 GeV, \l\ < 30) = J F(D)dN SNR 

-f 

Jd, 



' — I (HasMAlD) 



i=- 

1 6 



6.6 x 10 9 photons cm" 



2 s- 1 sr- 1 



(16) 



the source luminosity predicted for HESS sources would be 6.5 x 10 33 ergs/s and their 
average distance 3.5 kpc. Fig 2 shows the values of luminosity L and fraction of 7-ray 
emitting pulsars or SNRs, allowed by the condition that HESS detects 18 sources. 
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Fig. 2. — Values of luminosity L and fraction of 7-ray emitting pulsars or SNRs allowed by 
the condition that HESS detects 18 sources. The shaded region is not allowed. 



In the HESS catalogue of fifteen new sources 8 sources have SNRs as possible coun- 
terparts, and 4 of these new sources are potentially associated to pulsar wind nebulae. So 
pulsars and SNRS both occur but their surface densities are likely not so different. 



3. Contribution of unresolved sources to Milagro diffuse galactic emission 

Milagro is a water Cherenkov telescope, which has surveyed most of the northern sky 
at TeV energies (Atkins et al. 2000, 2003). Milagro has measured TeV energy fluxes from 
the Crab (Atkins et al. 2004) and the Cygnus Region (Smith et al. 2005) and, recently 
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Fig. 3. — Pulsar and SNR surface distributions versus the heliocentric distance D for the 
regions of the Galactic plane surveyed by Milagro ( 40° < / < 100° ) and by HESS ( -30° < 
/ < 30° ). 



the diffuse Galactic 7-ray emission up to TeV energies. Milagro Galactic diffuse emission 
is (7.3 ± 1.5 ± 2.3) x 1CT 11 photons s _1 cm' 2 sr _1 for E > 3.5 TeV in the region 
40° < / < 100°. (Atkins et al. 2005) The detection of TeV 7-ray diffuse emission is of 
fundamental importance in solving the long standing problem of the production of high 
energy 7 radiation in our Galaxy. In fact the Milagro result seems to exclude an additional 
hard spectrum component continuing to above 10 TeV to explain the 60 per cent excess of 
EGRET flux compared to 7r -decay production mechanism due to the local cosmic ray flux. 
(Hunter et al. 1997) In order to discriminate among different production mechanisms of the 
diffuse radiation, it is important to estimate the amount of the measured diffuse emission 
which could arise from unresolved pointlike sources. The Milagro sensitivity to point sources 
is about 10~ n photons s _1 cm -2 for E > ITeV and, in fact, it is of the order of 
the brightest integral source flux detected by HESS. However, part of the diffuse emission 
measured by Milagro is emission from unresolved HESS-like pointlike sources. 

The integral flux expected from the Milagro field of view from HESS-like sources as- 
suming a pulsar shaped surface density and source luminosity L = 1.7 x 10 33 erg/s is 



F P uisar{E > 200 GeV, 40 < / < 100) = J F(D) dN } 



pulsar 



bit 

!) 



30fcpc jL rl 

0.3 kpc 

2.4 x 10~ 9 photons s -1 cm" 2 sr _1 



dDD— J ^ dla pu i sars (l,D) 



9 
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(17) 

while from HESS-like sources assuming a SNR shaped surface density and L = 6.5 x 
10 33 erg/s is 



F snr (E > 200 GeV, 40 </< 100) = f F ( D ) 

= dDD — dla SNR (l,D) 

J 0.3 kpc U Jl=%- 



30 kpc rl= 
3 kpc 

3.1 x 10~ 9 photons s -1 cm^sr -1 , 



(18) 



The fluxes in Eq. (17) and (18) are lower than HESS flux from resolved and unresolved 
sources. In fact as shown in Fig.3 HESS observes the inner region of the Galaxy, whereas 
Milagro field of view is more spread toward the outer regions in the Galaxy where the 
number of sources is substantially lower, so a lower contribution from unresolved sources is 
expected for Milagro's region of the Galactic plane than for HESS. 

Milagro Galactic diffuse emission is measured for a threshold energy of about 3.5 TeV, 
whereas the fluxes from unresolved sources calculated in Eq.(17) and (18) refer to HESS 
threshold energy of 200 GeV. In order to estimate the contribution of unresolved sources to 
Milagro diffuse emission we correct the flux in Eq.(17) and (18) for the different threshold 
energy. All HESS sources are fitted by a power law spectrum 

<b{E,T) = %{-^) (19) 

and the average spectral index T = 2.32. The integral flux from unresolved pulsars if 5 per 
cent of radio loud pulsars emit 7-rays corrected for Milagro threshold becomes 

FpuisariE > 3.5 TeV, 40 < I < 100) = F pulsar (E > 200 GeV, 40 < I < 100) 



x 



F (-r+i) 

• -'-'Milagro \ 

Ehess 

= 5.5 x 10 -11 photons s~ x cm~ 2 sr _1 



If all HESS sources coincide with SNRs 

F snr (E > 3.5 TeV, 40 < I < 100) = F pulsar (E > 200 GeV, 40 < I < 100) 



(20) 
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x 



p (-r+i) 

- - C/ Milagro \ 



HESS 



7.1 X 10 



-11 



photons s 1 cm 2 sr 1 



(21) 



where Eunagro and Eh ess are Milagro and HESS energy threshold, respectively. The contri- 
bution from unresolved HESS-like sources to Milagro diffuse emission is comparable to the 
diffuse emission itself. Despite the uncertainty of this calculation, concerning especially the 
ratio bewteen 7-ray loud pulsars and SNRs, the result seems to indicate that unresolved 
point sources contribute significantly to Milagro diffuse 7-ray emission. Fig. 4 shows the 
fraction of the measured Milagro flux which unresolved sources contribute to for different Lo 
and for various total numbers of Galactic SNRs and pulsars. The lines in blue represent the 
conditions for HESS to detect 3, 9 and 18 sources, respectively. 
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Fig. 4. — The red lines represent the fraction of the measured Milagro flux which resolved and 
unresolved sources contribute to for different Lq and for different total number of Galactic 
SNRs and pulsars which are 7-ray loud. The lines in blue represent the conditions for HESS 
to detect 3, 9 and 18 sources, respectively. 



4. Contribution of HESS sources to EGRET diffuse galactic emission 

EGRET surveyed the 7-ray sky in the energy region between 30 MeV and 100 GeV. 
The Third EGRET Catalog lists 80 sources for latitudes less than ±10 degrees, of which 
5 are pulsars, one is a solar flare and 74 are unidentified sources. Studies of the spectral 
characteristics (Merck et al. 1996) and of the luminosity function (Mukherjee et al. 1995) 
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of EGRET unidentified sources show as unlikely that even a small part of these sources 
are pulsars or SNRs. The contribution of these sources to EGRET flux was evaluated 
by EGRET Collaboration and should be less than 10 per cent. Recently Strong (2006) 
showed that the contribution of unresolved sources with luminosity L(E > lOOMeV) < 
1 x 10 35 photons s -1 , which are too dim to be detected by EGRET, could provide up to 
5, 10 per cent of EGRET emission. Also he showed that combinations of different source 
populations could explain the excess 7-ray emission measured by COMPTEL and EGRET 
at MeV and GeV energies, respectively. 

The HESS new population of sources does not seem to coincide with EGRET Catalogue 
(Aharonian et al. 2006). From Gonthier et al. (2002) the EGRET point source sensitivity 
in the Galactic plane is about 1.6 x 10 -7 photons s _1 cm" 2 for E > 100 MeV. If 
we can believe that HESS TeV source fluxes continue down to GeV energies, ten of these 
sources would have been detectable by EGRET if extrapolated without any cut-off to 100 
MeV. The fact that only a few of HESS new sources are spatially consistent with the 95 per 
cent positional error of unidentified EGRET sources ( HESS J1640-465, HESS J1745-303 
and HESS J1825-137 for which a simple power law estrapolation of the spectrum to 1 GeV 
gives EGRET flux ) suggests that the simple estrapolation of HESS source fluxes to EGRET 
energies is probably not correct. These sources might have a cut-off close to 10 GeV. All 
these sources will be detectable by GLAST, which will have a flux threshold of about 5.0 x 
10~ 9 photons s _1 cm" 2 . The sum of the integral fluxes of the three HESS sources, which 
could coincide with EGRET unidentified sources, is 1.1 x 10~ 6 photons s _1 cm" 2 sr _1 
at 1.4 GeV, about one per cent of the EGRET diffuse flux. 

Following the method introduced in Chapter 2 we evaluate the contribution due to 
HESS-like sources to EGRET diffuse emission F sources (E > E , -30° < I < 30°) at different 
energies E above 0.50 GeV and for the two cases in which either all SNRs or 5 per cent of the 
radio loud pulsars contribute. In Table 1 the first column gives the energy E above which 
the integral flux F sources (E > Eq, —30° < I < 30°) is evaluated. The measured EGRET 
integral flux for E > Eq is given in the second column. The contribution of HESS-like 
pulsars or SNRs is given in columns third and fourth, respectively. As shown in Table 1 
the contribution from HESS-like sources can account for about 10 per cent of the GeV 
excess. This contribution becomes more important at higher energies. At about 85 GeV the 
contribution from either HESS-like pulsars or SNRs is about 20 per cent of the diffuse flux. 
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1.0 x 10~ 6 (8%) 
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2.3 x 10~ 6 


3.5 x 10- 7 (15%) 


3.3 x 10- 7 (14%) 


15.0 


1.4 x 10~ 6 


1.3 x 10- 7 (9%) 


1.2 x 10- 7 (9%) 
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4.4 x 10~ 7 


4.2 x 10- 8 (10%) 


4.0 x 10- 8 (9%) 
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6.4 x 10~ 8 


1.3 x 10~ 8 (20%) 


1.2 x 10~ 8 (19%) 



Fig. 5. — Table 1. Integral fluxes for EGRET measured diffuse 7-ray spectrum of the 
inner Galaxy region, (—30° < I < 30° and —5° < b < 5°) compared with the contribu- 
tion from HESS-like pulsars and SNRs for different energies. All fluxes are in units of in 
photons s _1 cm -2 sr _1 . EGRET fluxes are from Strong et al. (2005) . 



5. Conclusions 

Using HESS new sample of Galactic 7-ray emitters we estimated the contribution of 
unresolved sources to the diffuse emission measured by EGRET and by Milagro. We find 
that the contribution of unresolved HESS-like sources to the diffuse 7 ray emission is large at 
energies above 100 GeV, whereas it amounts to about 10 per cent of EGRET measurements 
of the diffuse flux at GeV energies. Since pulsars and SNRs are distributed mostly in the 
Galactic plane, their unresolved fluxes will mainly contribute to the low latitude diffuse 
emission. Thus we expect that the total flux due to unresolved sources contributes less at 
higher latitudes. 

The main uncertainties of this calculation consist in assuming that all sources are mono- 
luminous and the distribution of 7-ray sources follows the distribution of either pulsars or 
SNRs. Of the 1509 pulsars catalogued by the Australian National Telescope Facility (Manch- 
ester et al. 2005) only a few are also 7-ray emitters. This could be due to strong selection 
effects as suggested by the existence of pulsars, such as Geminga, which are 7-ray loud pul- 
sars, yet not observed in the radio. However, both the SNR and pulsar distribution give 
similar conclusions, as would any class of sources that follows a similar spatial distribution. 

To draw more certain conclusions about the high energy 7-ray Galactic sky, new obser- 
vations are of fundamental importance. VERITAS will be able to survey the same region 
of the sky as Milagro. Milagro's recently improved sensitivity will better image the dif- 
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fuse emission. Finally, GLAST will investigate the window of energy between 10 MeV to 
300 GeV, covering the energy gap left between EGRET and the ground based gamma-ray 
observatories. 

We would like to thank Andrew Strong, Igor Moskalenko and Nino Panagia for useful 
discussions and comments. 
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